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1. Introduction 

The first-order phase transition plays an important role in diverse fields of physics. Its 
most characteristic feature is that a phase boundary is a coexisting region, where two or more 
phases coexist. The dimension D of a coexisting region in the phase diagram depends on the 
number r of coexisting phases. 1-5 It also depends on which variables are taken as the axes of 
the phase diagram. 1-3 

For example, for a single-component system whose natural variables of the entropy 1-4 are 
U (energy), V (volume) and N (amount of substance), the phase diagram may be drawn 
either in the T-P plane, or in the T-v plane, or in the u-v plane, where T, P, v and u are 
the temperature, pressure, molar volume (= V/N) and molar energy (= U/N), respectively. 
Here, v and/or u are used instead of V and/or U in order to make the diagrams independent 
of N. These phase diagrams are schematically shown in Figs. 1, 2 and 3, respectively. 3 We 
denote D for these diagrams by D(T,P), D(T,v) and D(u,v), respectively. For example, 
D(T,P) = l,D(T,v) = D(u,v) = 2 for a liquid-gas coexisting region, where two phases 
coexist (i.e., r = 2), whereas D(T,P) = 0,D(T,v) = l,D(u,v) = 2 at the triple point, where 
r = 3. 1-3 Since T and P can take constant values across a coexisting region, 1-3 D of the 
coexisting region tends to shrink if T and/or P is taken as an axis(es) of the phase diagram, 
and thus D(T,P) < D(T,v) < D(u,v). 

For such a single-component system, D(T, P) coincides with the 'thermodynamical degrees 
of freedom' /, which is given by the Gibbs phase rule 1-5 (see eq. (7) below). On the other hand, 
D(u,v) > D(T,v) > f in general. 1 ' 3 For example, / = 1 = D(T,P) < D(T,v) = D(u,v) = 2 
for a liquid-gas coexisting region, whereas / = = D(T,P) < D(T,v) < D(u,v) = 2 at the 
triple point. 

The situation becomes more complicated for a multi-component system, which con- 
sists of q (> 2) different substances. We assume that the natural variables of entropy are 
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Fig. 1. A schematic phase diagram of a single-component system, plotted in the space spanned by the 
axes corresponding to T and P. The solid- liquid-gas coexisting region is a point, which is called 
the triple point. On the other hand, the liquid-gas, solid-liquid and solid-gas coexisting regions 
are lines, which are called coexisting lines. 



U,V, Ni, . . . , Nq, where Nj (j = 1, . . . , q) is the amount of substance j. That is, the entropy 
S is a function of these additive (extensive) variables; 

S = S(U,V,N u ...,N q ), (1) 

which is called the fundamental relation. 1, 3,4 Let A^ot be the total amount of substances, 

N tot = N 1 + ... + N q . (2) 



To make phase diagrams independent of iVtot, we use the normalized variables; 

U 

V 



tot 



(3) 
(4) 



N J- (J = l,. •.,<?). (5) 



iVtot 

Among m, . . • , n q , only (q — 1) variables are independent because of the identity 

2>j = l. (6) 

i=i 

Hence, the phase diagram may be drawn in the (g+l)-dimensional space that is spanned by the 
axes corresponding to either (T, P, m, . . . , n q -i), (T, v, nx, . . . , n q -i), or (u, v, m, . . . , 
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Fig. 2. A schematic phase diagram of a single-component system, plotted in the space spanned by 
the axes corresponding to T and v (= V/N). Unlike Fig. 1, the solid-liquid-gas coexisting region 
(triple point) is a vertical line, whereas the liquid-gas, solid-liquid and solid-gas coexisting regions 
are planes. 



which correspond to the so-called TPN, TVN and UVN representations, respectively. 
We denote the dimensions of a coexisting region in these spaces by D(T, P, m, . . . , n q -i), 
D(T, v, ni, ■ ■ ■ , fiq-i) and D(u, v, m, . . . , n g _i), respectively. 

On the other hand, the thermodynamical degrees of freedom 1-5 / is defined as the number 
of variables that can be varied independently in the coexisting region, among the intensive 
variables T,P, fix, . . . ,fi q , where (ij denotes the chemical potential of substance j. The Gibbs 
phase rule gives 1-5 

/ = q + 2 - r. (7) 

For q > 2 neither of -D's coincides with / in general, and the explicit formulas for D's were 
unknown. 

Knowledge about D's would be very helpful in drawing phase diagrams of new materi- 
als. 1-5 Since phase diagrams are most fundamental to studies of macroscopic systems, gen- 
eral formulas would be valuable which give D(T, P,m, . . . , n q -i), D(T, v, ni, . . . , n q -i) and 
D(u, v, ni, . . . , n q -i) as functions of q and r. The purpose of the present paper is to derive 
such general formulas. We also evaluate D's in some other spaces. The results will be sum- 
marized as the chain of equalities and inequalities, eq. (41). 
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solid+liquid+gas (triple point) 



V 



u 



Fig. 3. A schematic phase diagram of a single-component system, plotted in the space spanned by 
the axes corresponding to u (= U/N) and v (= V/N). All the coexisting regions are planes. For 
example, the triple point, which was a point in Fig. 1, extends to the shaded area in this figure. 
Unlike Figs. 1 and 2, every point in this diagram corresponds to a single equilibrium state, and 
different points correspond to different equilibrium states. 

2. UVN representation 

2.1 Basic relations 

We consider a g-component system, the natural variables of whose entropy 1-4 are assumed 
to be U,V, Ni, . . . , N g . For each coexisting region, we label the coexisting phases by a = 
1, . . . , r. For example, the shaded region of Fig. 3 is a coexisting region, in which the gas 
(a = 1), liquid (a = 2) and solid (a = 3) phases coexist, hence is called the triple point. Note 
that there is the trivial r!-fold degeneracy in labeling a = 1, . . . , r. Since this degeneracy does 
not affect the values of -D's, we henceforth forget about it. 

The values of U, V and Nj in phase a are denoted by U a , V a and N", respectively. They 
are related to those of the total system as 



r 



u 




(8) 



a=l 



r 



V 




(9) 



a=l 



r 




(10) 



a=l 
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Corresponding to u, v, nj of the total system, we define 



jja 

Jv tot 
ya 

J\f a ' 



u = 

a=l 



a=l 



(11) 

(12) 



n« - (j = l,...,q) (13) 



for each phase, where N^ t is the total amount of substances in phase a, 

N t a ot ^N 1 * + ... + N q * (a = l,...,r). (14) 
Corresponding to eq. (6), we have 

q 

J^raJ = 1 for all a. (15) 

3=1 

We also define the molar fraction of each phase 

N a 

^ = jf- (a = l,..-,r) (16) 

which satisfies the trivial identity 

E^ = l- (17) 

a=l 

By dividing eqs. (8)-(10) by -/V to t, we obtain 

r 

J2" a u a , (18) 



v = ^i/V, (19) 



£>X (i = !>•••><?)• (20) 



nj = 

a=l 

Note that rac/i phase is, by definition, 3 ' 5 ' 6 homogeneous spatially. Hence, the local value of 
the molar energy is equal to u a at everywhere in phase a. On the other hand, the total system 
is inhomogeneous in a coexisting region, and the local value of the molar energy is equal not 
to u but to either of u 1 , u 2 , . . ., or u r . As seen from eq. (18), u actually represents the average 
molar energy, which is the weighted average of u a, s of the coexisting phases. The same can 
also be said about v and rij. 

Intensive variables in phase a are continuous functions of u a , v a , nf, . . . , n"_ 1 , 1 ' 3 and have 
the same values as those of the total system; 1-4 

T = T(u a , v a ,n?, = T(u, v,m,..., n,_i), (21) 

P = P(u a , v a ,nf, nj_!) = P(u, v,m,..., n q ^), (22) 

Hj = fij(u a ,v a ,nf,. . . ,71^) = n j {u,v,n 1 ,...,n q - 1 ) (j = l,...,q). (23) 
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Here, the functions T(-, •,...,•), P(-, •,...,•) and fij(-, -,...,-) are obtained by differentiating 
the fundamental relation, eq. (I). 1 ' 3 ' 4 

2.2 D in the {u a , v a , nf , . . . , n^_ 1 , ^ a } a =i,...,r space 

Our principal purpose is to evaluate D's in the (q + l)-dimensional spaces such as the one 
spanned by the axes corresponding to (T, P,m,..., n q -\). As a preliminary, we first evaluate 
the dimension D({u a ,v a ,n" , . . . , n° 1 , v a } a ) in a larger space which is spanned by the axes 
corresponding to r(q + 2) variables {u a ,v a , nf, . . . , n q l _ l , v^ ot ] a= i ^ r . 7 

Since r different phases coexist, we have 

T{u\v\n\, n^) = . . . = T(u r ,v r ,n r l , n^), (24) 

P{u\ v 1 ,n{,..., nj_!) = . . . = P(n r , , . . . , n r q _ ± ), (25) 

^•(u 1 . . . .n^i) = . . . = ^j{u r ,v r ,n\, . . .,71^) (j = l,...,q), (26) 

which impose (r — 1) x (<7 + 2) conditions on r x (q+ 1) variables {u a , v a ,n", . . . , rag_ 1 } Q= i ) ... ir . 
Therefore, the dimension of the set of values of {u a ,v a ,nf, . . . , n^_ 1 } Q= i ) ... )r is evaluated as 

r(q+l)-(r-l)(q + 2)=q + 2-r. (27) 

Among the residual variables is 1 , . . . , u r , we have eq. (17). Hence, 

D({u a ,v a ,n?,...,n°_ 1 ,v a } a ) = (q + 2 - r) + (r - 1) 

= g + (28) 

2.3 -D in i/ie u,v,ni, . . . , n q _\ space 

With the help of eq. (28), we can evaluate D(u, v, m, . . . , n q -\) as follows. It is ob- 
vious from eqs. (8)-(10) that the values of U, V, Ni, . . . , N q are uniquely determined by 
the values of {U a , V a , , . . . , iV°} a= i v .. jr . Furthermore, the latter (including the num- 
ber r) is uniquely determined by the former, because otherwise two different states would 
have the same value of the total entropy S(U, V, Ni, . . . , N q ), and thus the total sys- 
tem would be unstable. 3 Therefore, the values of {U a , V a , , . . . , iV°} a= i v .. jf . (including 
r) have one-to-one correspondence with the values of U, V, Ni, . . . , N q of the total sys- 
tem. This means that {u a , v a , n", . . . , A r t ot}o=i, nave one-to-one correspondence with 
u, v, ni, . . . , n q -i, A?tot- Hence, {u a , v a , . . . , n^_ 1: v a } a =i,...,r have one-to-one correspon- 
dence with u, v, n±, . . . , n q -\. Therefore, 

D(u,v,n u ----,n q -i) = D({u a , v a , n q _ 1: u a } a ) 

= 9 + 1, (29) 

which is equal to the dimension of all possible states for a given value of AW 1 - 3 - 4 This is rea- 
sonable (obvious in some sense) because each equilibrium state has one-to-one correspondence 
with the set of values of U, V, Ni , . . . , N q , even when the equilibrium state is in a coexisting 
region. 1,3 
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For a single-component system, for example, formula (29) gives D(u,v) = 2, which is 
consistent with Fig. 3. 

3. TVN representation 

3.1 D in the T, {v a ,n®, . . . , n^_ 1 } a= i ! ... ;r space 

To evaluate D(T,v,ni, . . . ,n g _i), we start with considering D in a larger space which is 
spanned by the axes corresponding to (T, {v a , nf, . . . , n^_ 1 } Q= i ) ... jr ). 7 

Since r different phases coexist, we have 

P 1 (T, v\ n\, . . . , «£_!) = . . . = P r (T, v r , n{, . . . , n*^), (30) 
t x](T,v\n 1 1 ,...,n 1 q _ 1 ) = ... = » r j (T,v r ,n r 1 ,...,n r q _ 1 ) (j = l,...,q), (31) 
which impose (r — 1) x (q + 1) conditions on (1 + r x q) variables. Therefore, 
D(T,{v a , = l + rg -( r -l)(g + l) 

= g + 2-r, (32) 

which coincides with /. 

D in i/ie T, w, m, . . . , n g _i space 
With the help of eq. (32), we can evaluate D(T, v,ni,..., n q -\) as follows. Equations (19) 
and (20) show that for each set of values of (T, {v a , nf , . . . , n"_ 1 } Q= i ) ... )T .) one can vary q 
variables v,ni, . . . , ra ? -i by varying r variables {^ a } a =i,..., r subject to one condition, eq. (17). 
Therefore, 

D(T,v, ni ,...,n q -i) = D{T,{v a ,n'},...,n q t _ 1 } a )+mm{q,r-l} 

= q + 2 - r + min{g,r - 1} (33) 

'«+! (r< ?+ D, m 
q (r = q + 2), 

where in the last line we have taken account of an important consequence of eq. (7); 

r < q + 2, (35) 

which gives the upper limit of the number of coexisting phases. 1-4 

For a single-component system (q = 1), for example, the above formula gives 

D(T,v) = < 



2 (r<2), 
1 (r = 3), 



(36) 



which is consistent with Fig. 2. 
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4. TPN representation 

4-1 D in the T, P, {n", . . . , n? : _ 1 } Q;= i ) ... )r space 

To evaluate D(T, P, ni, . . . , n q -±), we first consider D in a larger space which is spanned 
by the axes corresponding to (T, P, {n", . . . , rag_i}a=i,..., r )- 

Since r different phases coexist, we have 

fi 1 j (T,P,n 1 1 ,...,n 1 q _ 1 ) = ... = f i r j (T,P,n r 1 ,...,n q _ 1 ) (j = l,...,q) (37) 

which impose (r — 1) x q conditions on [2 + r x (q — 1)] variables. Therefore, 

D(T,P,{nf,...,nJ_ 1 } a ) = 2 + r(g - 1) - (r - l)g 

= g + 2 - r, (38) 

which coincides with /. This well-known result is also called the phase rule. 3 ' 4 
4-2 D in the T, P,n±, . . . , n q _i space 

It is seen from (20) that for each set of values of (T, P, {n®, . . . , rag_i} a =i,..., r )) on e can 
vary (q — 1) variables ni, . . . , n g _i by varying r variables {v a } a =\,...,r subject to one condition, 
eq. (17). Therefore, 

D(r,P,ni,...,ng_i) = J D(T,P,K,...,n°_ 1 } a )+min{g-l,r-l} 

= g + 2-r + min{g-l,r-l} (39) 

q+l (r<q), 
= < V _ ' (40) 

2g+l-r (g + l<r<g + 2), 
where in the last line we have taken account of eq. (35). 

For a single-component system (q = 1), for example, this formula gives D(T,P) = 3 — r, 
which coincides with / and is consistent with Fig. 1. 

5. Conclusions 

Our principal results are eqs. (29), (34) and (40). We have also derived additional results, 
eqs. (28) and (32). For completeness, we have also described the known results, eqs. (7) 
and (38). Here, eq. (7) can also be written as D(T, P, fj,i, . . . , fi q ) = q + 2 — r, because, by 
definition, / is the the dimension of a coexisting region in the space that is spanned by the 
axes corresponding to (T, P, fii, . . . , fi q ). 7 

By collecting all these results, we obtain the following chain of equalities and inequalities; 

< D(T,P,fi 1 ,..., liq ) = q + 2-r = f 

= D(T, P, K, . . . , n£_J Q ) = D(T, K, nf , . . . , <_J a ) 

< D(T, P, m, . . . , n q -i) = q + 1 - r + minjg, r} 

< P(T, u, m, . . . , n q -i) = q + 2 - r + min{g, r - 1} 

< £>(«,«, ni,...,n,_i) =P(K, V a , <_ 1; ^} Q ) =q+l. (41) 
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Here, q + 1 is the dimension of all possible states for a given value of iVtot- 1 ' 3 ' 4 It also 
agrees with the dimension of the space that is spanned by the axes corresponding to ei- 
ther (T, P, n±, . . . , n q -i), or (T, v, ni, . . . , n q -i), or (u,v,ni, . . . , n q -\). This chain of equalities 
and inequalities may be regarded as the fundamental phase rule, which shows clearly how 
the dimension of a coexisting region varies depending on the choice of the variables. It will 
be helpful in studying first-order phase transitions and in drawing phase diagrams of new 
materials. 

Finally, we note the following points. Although we have assumed that the natural variables 
of entropy 1-4 are U, V, N±, . . . , N q , generalization to other cases (such as the case where they 
include the total magnetization 3 ' 4 ) is straightforward. Furthermore, we have assumed, as in 
the case of the Gibbs phase rule, that there is no accidental degeneracy among equations 
which have been used in calculating D's. Hence, it is in principle possible (though would be 
rare) that D's take values that are different from our formulas. 
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